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A method for studying initial geometry fluctuations via event plane correlations in 

heavy ion collisions 
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A method is proposed to measure the relative azimuthal angle distributions involving two or 
more event planes of different order in heavy ion collisions using a Fourier analysis technique. The 
analysis procedure is demonstrated for correlations involving two and three event planes ($ n , 5>m and 
$h). The Fourier coefficients of these distributions are found to coincide with previously proposed 
correlators, such as cos(6$2 — 6$3) and cos($i +2$2 — 3$3) etc, hence the method provides a natural 
framework for studying these correlators at the same time. Using a Monte Carlo Glauber model to 
simulate Au+Au collisions with fluctuating initial geometry, we are able to identify several new two 
or three-plane correlators that have sizable magnitudes and should be measured experimentally. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION 

In heavy ion collisions at the Rclativistic Heavy Ion 
Collider (RHIC) and the Large Hadron Collider (LHC), 
the fluctuations of the positions of nuclcons in the overlap 
region are found to play an important role in controlling 
the shape of the initial geometry of the created matter, 
which subsequently controls the azimuthal anisotropy of 
the particles in the final state [1-3]. The shape of the 
geometry in azimuth can be characterized by a set of 
multi-pole components (also known as "eccentricities") 
at different angular scale, calculated from the partici- 
pants and binary collisions at (r, <j)) [4] : 



\/ (r 2 cosncf)) 2 + (r 2 smruf)) 2 

V) ; 



(1) 



21 with a weight of S = 0.14 for binary collisions and 

22 (1 — S)/2 = 0.43 for participants [5], where (r,4>) are 

23 calculated relative to the weighted center of gravity. The 

24 orientations of the minor and major axes for each mo- 

25 ment n are given by 



$n 



atan2( (r 2 sin n<t>) , (r 2 cos n<j)) ) n 



and 



« = ®r 



■k jn 



(2) 



(3) 



27 respectively. The minor axis direction $„ is also known 

28 as the n -order participant plane (PP). The values of e n 

29 and <f> n can be calculated easily using simple geometric 

30 models such as Monte Carlo Glauber code from [6]. 

31 When fluctuations are small and linearized hydrody- 

32 namics is applicable, each e„ is expected to indepen- 

33 dently drive the corresponding n th -order anisotropic flow 

34 v n along $„ [4]. In this case, one may rely on a simple 

35 Glauber model calculation to estimate the correlations 



36 between anisotropic flows of different order . Previ- 

37 ous studies [3, 7-11] show that significant correlations 

38 can exist between $2 and $4 due to the almond shape 

39 of the average collision geometry. However, correlations 

40 involving odd planes for n > 2 are found to be gener- 

41 ally weak except in very peripheral collisions, e.g. be- 

42 tween $ 2 and $3 or between $ 2 an d <f> 5 [7]. Experi- 

43 mental studies support a strong correlation between <J> 2 

44 and $4 [12, 13], but a weak correlation between $2 and 

45 $3 [14, 15]. The correlations among three planes of diffcr- 

46 ent order have also been investigated recently [3, 10, 11], 

47 such as $1 + 2$ 2 - 3$ 3 and $! - 4$ 2 + 3$ 3 ; they are 

48 argued to contain strong correlations between the dipole 

49 asymmetry and the triangularity. Here we propose an al- 
so ternative experimental method for measuring these cor- 

51 relations. The expected performance of this method is 

52 evaluated based on the correlation signals from Glauber 

53 model. 



54 II. METHOD 

55 The n th -order flow has n-fold symmetry in azimuth, 

56 and the correlations between flow directions $„ and $ m 

57 are completely described by the differential distribution 

58 dN evts /(d (fc($n — $m)))) with k being the least common 

59 multiple of n and m, i.e. k = LCM(n, m). This distri- 

60 bution should be an even function due to symmetry, and 
6i can be expanded into a Fourier series: 



dN c - 



d(fc($„-$ m )) 

V 3 



l + 2j2vl m cosjk($ n -$ m l4) 



(cos jfc($„ - $ m )} 



(5) 



In a real experiment, the underlying true event plane 
directions 4> ra and <& m are unattainable due to limited 



1 This was found, via a hydrodynamics calculation, to be approx- 
imately true for v n < 3, but not so for n > 3 except in central 
collisions [9, 16]. Our estimation of higher order <J? n should be 
digested with this caveat in mind. 



detector acceptance and finite multiplicity. They are ap- 
proximated by the measured event plane angle *5/ n and 
$f m , calculated based on the azimuthal distribution of 
particles in the detector acceptance. The coefficients 
(cosjfc(<£> n — $,„)) can be obtained by calculating the 
raw coefficients (cos jk(^ n — ^ m )), followed by a simple 
correction for finite event plane resolution: 



3 are required (see Fig. 5). The coefficients are: 



V j (cosjfc(E n -tt m )) 

n ' m Res{jfc*„}Rcs{jfc* m } 
Res{jfc*„} = (cos jfc(*« - *«)) 



71 The resolution factor Res{jfc^ n } can be determined us- 

72 ing the standard two-subevent or three-subevent meth- 

73 ods [17]. To avoid auto-correlations, the $„ and $ m 

74 should be measured using sub-events covering different r\ 

75 ranges, preferably with a gap in between. 

76 Interestingly, some of the two plane correlators defined 

77 by Eq. 5 are related to the so called mixed harmonics, 

78 referring to vi n measured in <f>„ event plane for integter 

79 I > 2, denoted as vi n {& n } (see Ref. [17]). For example, 
so it is straightforward to show that V r 2 1 4 i s simply the ratio 
si of the integral v 4 measured in the <J> 2 plane (i^l*^}) 

82 to the integral v 4 measured in the $ 4 plane (tj 4 {$ 4 }): 

83 (cos4($ 2 — ^4)) = Vi{^2}/vi{^4}- More generally, one 

84 has: 



(cosm($„-$ m )) 



i.{$m} 



. mod n = 0. (8) 



K l,2.3 



(cos (z$ 2 ,i ± j$3,i)) 

(cos«$2,i cos j$3 ; i) =F (sinz$2,isinj$3 i i) (11) 



Under this notation, the two-plane correlator can be 



104 treated as special case: V x 



i,0 
2,3 



vu,vi 



0,3 _ 
2,3 



y/ 3 and 



105 

(6) 106 
107 

(7) 108 



V3 2 - The average of the sine term in Eq. 11 



T/3.7-2.7 
r l,2,3 

may not be zero, if the fluctuations of $2 and $3 relative 
to $1 are correlated, i.e. $2 and $3 prefer to appear si- 
multaneously to one side of $1 . It represents a non-trivial 
correlation that is of great interest for understanding the 
nature of the fluctuations (see also [3]). 

A similar resolution correction procedure can be used 
to connect the measured correlated with the corrected 



V l ' ±3 - 

^1,2,3 — 



(C0S(J*2,1± j*3,l)> 



Res{|2i ± 3j|*i}Res{2i* 2 }Res{3j* 3 } 



(12) 



as Additional examples include (cos6(<I > 2 — ^o)) = 
as v 6 {$2}M{$6} and (cos6($ 3 - $ 6 )) = v 6 {$ 3 }/v 6 {<£>(i}- 

87 The method described above can be generalized to 

88 correlations involving three or more event planes. As 
as pointed out in Ref. [10] , the correlations that can be mea- 

90 sured experimentally, involve combination of I planes of 

91 different order: ci$i + 2c2 ( I > 2--- + Ici^i with c\ + 2c2... + 

92 lei = 0. The correlations involving three planes of diffcr- 

93 ent order, e.g. $1, $2 and $3, have the following form: 



where we have assumed that \]/ n is distributed randomly 
around <f>„, such that (sin jn (^>„ — $n)} = 0. Again, the 
^1, ^2 and ^3 should be calculated from subevents cov- 
ering different 77 acceptances to avoid auto-correlations. 

In [3], Teaney and Yan proposed to study the correla- 
tor cos($i + 2$ 2 - 3$ 3 ) and cos($i - 4$ 2 + 3<f> 3 ). In 
our notations, they correspond to the cosine average of 
the 2-D distribution ($2,11^3,1) projected along the di- 
rection (i,j)=(l,-l) and (2,-1), respectively. Our frame- 
work provides a natural way to visualize and systematize 
the study of these correlators. 

Other triple plane correlators can be similarly ana- 
lyzed, the first few are 



Ci$i +2c 2 $ 2 + 4c 4 <I>4 
c 1 <& 1 + 3c 3 $ 3 + 4c 4 $4 

202*2 + 3c 3 $3 + 4c 4 $ 4 

ci$i +2c 2 $2 + 5c 5 $ 5 

Ci$! +3C2$3 + 5C 5 $ 5 



C2$2,l+C 4 $4,l 

C 3 $3,l+C 4 $4,l 
C3 

2 

C2$2,l 



3,2 



(13) 

(14) 

C 4 $4, 2 (15) 
C5$5,l (16) 



C3$3,l + C 5 $5,l (17) 



C22($2-*l)+C 3 3($3-$l), 
C2$2,l +C 3 $3,1 (9) 



ci$i +2c 2 $2 + 3c 3 $3 

where we use the constraint c\ + 2c2 + 3c3 = and 



95 we adopt the short-hand notations: $ ? 

96 <J> m ),* 



k($ r 



m J ? ^ 71, m 



k(^ n — 'I'm). This type of correlations 

97 can be generally determined from the underlying 2-D dis- 

98 tribution in ($2,1)^3,1) via a similar Fourier expansion 

99 approach: 



d$ 2 ,id*3,i 



oc 1 + 2 



OO 

Yl N'.a cos J $ 2,i + Vl 3 cos j$ 3 ,i 



i=i 



2 Yl ^3 COS (^2,l±i$3,l). (io) 



100 This series is expected to converge quickly for non- 143 

101 peripheral collisions. Therefore, only the terms for i,j< 144 



Note that C3/2 in Eq. 15 is an integer by the require- 
ment 2c2 + 3c3 + 4c 4 = 0. These correlators can be 
uniquely identified with one of the Fourier coefficients 
in the double differential distributions similar to Eq. 10. 
However, the expression of triple plane correlator in 
terms of the correlation between two di-plane corre- 
lators is not always possible, for example (cos(2$ 2 + 
3$3 — 5$5)). In this case, it can be regarded as a 
sum of the Fourier coefficients for triple differential dis- 
tributions (i 2 iV cvts /((i$3 ! 2d c E > 5.irf ( I > 5,3) that contribute to 
cos(2$ 2 + 3$3 -5<J> 5 ). 

The measurement of correlations involving two or more 
event planes are feasible at the LHC due to the large de- 
tector coverage in 77 (-5 < 7*7 < 5 in ATLAS and CMS), 
and excellent reaction plane resolution [18, 19]. This al- 
lows the choice of many non-overlapping sub-events, each 
with very good 77 coverage for these multi-plane correla- 
tion measurements. 



!■!'., 



The coefficients 7„ 3 m or K 



i,±j 



n ,m 



n.m,< 



146 previousiy proposed multi-particle correlators from 

147 Ref. [10, 11]. That approach effectively applies 

us a |c„|-particle weight 

149 if the angle nc n <& n appears 



are related to the mn event plane measurement, while the calibration proce- 
dure could be more involved for multi-particle correla- 
tions [23]. 



150 This weight 

151 and reduce 

152 have small 

153 cos(c„n$„ 4 
J 

154 V n 



1} {|C 



1} 



maximizes the 
the contribution 
v n . For cos(c„n<f>„ 

c m m$ m + c h h$ h ), the 



«A = (V„)l(v n )2-{v n )\c n \ 

in the correlator. 

correlation signal 

for events which 

m$ m ) and 

are 



and vi M} vJt ml} v { h {chl \ 



they 



1} {\Cr. 



then divided 

iyi^i>\ 



by 



,{l 



.iyi 



weights 
respectively; 
} ) and 

v h "") t° obtain the true correlations. 

157 Note that the weighting procedure can also amplify con- 

158 tributions from the tail of the e„ distribution, especially 

159 for large values of c„, this may complicate the mapping 
i6o from the measurement to correlations between e„. 

i6i In contrast, all events have the same weight in our 

162 approach (including those with small e„ values unfortu- 

163 natcly). In our opinion, the two methods arc compli- 

164 mentary to each other. In fact, it is possible to construct 

165 some hybrid correlators that involves azimuthal angle of 

166 both event planes and single particles. For example, one 

167 can consider the following mixed correlator between a + b 

168 particles and event planes \&„, \f m : 



(cos(nc n $ n — mCm $ m )} {q} {b} = 

(cos (J2 K' b m + n{c n - a)* n + m(c m - 6)* m )> 
{vi a } «m } }Res{n(c n - a)* rl }Res{m(c m - 6)* m } 

^2 ^™." 1 = n ( ( t >1 + ■• + ^ a ) + m {<t>a+i + ■■ + 4>a+b) (19) 



(18) 



(OOS6($2 - $3)) {2} 

(cos2($i-$ 2 )} (lOT)i){2} 
(cos($i+2$ 2 -3$ 3 ))«n>i 



(cos (30i +302 -6» 2 )) 
((«3)i(« 3 )2)Res{6* 2 } 

(COS (01 +02 -2^2)) 

((it)Wi)i(TODi) 2 >Res{2*2} 

(cos (0i +21/2 -2,^z)) 

(u.i;i)Res{2*2}Res{3*3} 



(20) 



(21) 



(22) 



175 where w(pt,ij) = px — (p|>) (v) I (pt) (v) ^ s the px and rj 

176 dependent weight (w is rapidity-even for Au+Au colli- 

177 sions) that designed to suppress global momentum con- 

178 servation effects and maximizing the v\ signal (or effec- 

179 tively increasing the resolution for ^1) [10, 15, 20, 21]. 

180 Even though terms related to v\ {<j}\ and/or fa) appear 
i8i in Eq. 21-22, the global momentum conservation effects 

182 are expected to be either negligible (Eq. 21) or absent 

183 (Eq. 22) for these correlators [15, 22]. 

184 The hybrid correlators are useful in real experiments 

185 where detector subsystems have limited geometrical ac- 

186 ceptance, finite granularity (so can't distinguish individ- 

187 ual particles), limited px reach or no px information at 
las all. In this case, it is straightforward to calibrate the 



192 

193 



III. EXPECTED BEHAVIOR BASED ON A 
GLAUBER MODEL 

A. Correlation between two planes 



1 I 0-5% j 



4(* 2 *-* 4 *) 




4(* 2 *-V) 



■ 0-5% 
45-50% 



169 where nci — mc rn = 0, fa, ...,0 a+ b are azimuthal angles 

170 of a + b particles, and subscript v„ »i indicates the 

171 weighting factor introduced by those particle multiplcts. 

172 Similar formula can be generalized to correlations involv- 

173 ing more than two event planes. Three interesting exam- 

174 pies are: 



FIG. 1: (Color online). (Top panels) The distribution of the 
angle difference between major axes for £2 and £4 in two cen- 
trality intervals. (Bottom panel) The Fourier spectra. 

195 Similar to [7], we use a simple Glauber model to esti- 

196 mate the level of the correlations between n th - and m th - 

197 order participant plane. About 2.5 Million Au+Au col- 

198 lisions are simulated, where each Au ion is populated 

199 randomly with nucleons with a hard-core of 0.3 fm in 

200 radii, according to the Woods-Saxon distribution with a 

201 radius of 6.38 fm and diffuseness of 0.535 fm. A nucleon- 

202 nucleon cross-section of a — 42 mb is used. However, 

203 instead of using the minor axes of e„ as the proxy for the 

204 true event planes, we actually calculate the correlations 

205 between the major axes, which are related to the former 

206 by a simple phase shift 8 n>m : 

HK-^n) = H$n ~ *m) + 5 n , m (23) 

5 n ,m = k(l/n-l/m)7r (24) 



207 The corresponding Fourier coefficients are denoted as 
cosfc($* - $£J), and are related to V^ m as 



yj* 



vr m = (cos sy vi, 



(25) 



209 The reason for doing this is a simple matter of conve- 

210 nience: the correlations between major axes are found to 

211 be almost always positive in the Glauber model, hence 

212 using major axes simplifies the presentation. It is intcr- 

213 esting to note that the phase shift, when folded to [0, 2tt], 

214 is (Sn tm mod 2tt) = or it. The latter case leads to a sign 



flip: 
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-V 3 . 
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FIG. 2: (Color online) Centrality dependence of first Fourier 
coefficient of the correlation V n * m = (cos fc(\P^ — ^m)) f° r 
various choices of n and m. 



216 The top panels of Fig. 1 show the 4($2 ~~ ^1) cor- 

217 relations predicted by the Monte Carlo Glauber model. 

218 The correlation is weak in central collisions but is quite 

219 strong in peripheral collisions. This is understood [7, 9] 

220 due to a detailed interplay between the fluctuation and 

221 average shape for the collision geometry: the central col- 

222 lisions are fluctuation-dominated, hence $„ are largely 

223 uncorrelated, while the peripheral collisions arc domi- 

224 nated by the average geometry which has tin compo- 

225 nents that are aligned [24]. Figure 1 also shows that 

226 the first order component captures most of the shape 

227 information in central collisions. In contrast, many com- 

228 ponents are needed to describe the tight correlation in 

229 peripheral collisions. This behavior is generally true in 

230 the Glauber model: whenever the first term V£* m is large, 

231 more higher-order terms are needed to describe the full 

232 distribution. Note that if the participant planes are used 

233 instead, the distribution 4($2 — $4) would show an anti- 

234 correlation: the distributions have their minima at 0, and 

235 the sign of V% 4 alternates between positive and negative: 

236 vi A = (-lyvfr. 

237 The calculations are extended for all types of corrcla- 

238 tions for k up to 16. Of course, additional correlations 

239 can also be calculated but the resolution is expected to 

240 deteriorate quickly for large values of k. The centrality 

241 dependence of these correlations, characterized by the 

242 first Fourier coefficient V^* m , are shown in Fig. 2, where 

243 the centrality is characterized by number of participat- 

244 ing nucleons (A^ part ). In most cases, the correlations are 

245 found to be either consistent with zero or positive (except 

246 for Kfg in mid-central collisions). In particular, strong 



247 correlations are observed for 4($2 — $4), 6(^3 - <&g) and 

248 6 ($2 — ^6 )i they are presumably associated with the av- 

249 crage geometry. The correlations are small in central col- 

250 lisions and over the full range for other choices of n and 

251 771, suggesting that the correlations are generally weak 

252 when they are dominated by fluctuations. Wc would 

253 like to draw the reader's attention to the bottom pan- 

254 els, which suggest that there are significant correlations 

255 between $1 and all other higher-order PPs. Recently, sig- 

256 nificant dipolar flow V\ has been observed in Pb-Pb col- 

257 lisions at the LHC by the ATLAS Collaboration [18] and 

258 a theoretical group [25] based on the ALICE data [26]; 

259 large dipolar flow is also predicted in hydrodynamic [27] 

260 or transport models [21]. Therefore, it is reasonable to 

261 assume that the correlations between dipolar flow and 

262 higher-order flow is large and measurable, as long as one 

263 can find a way to determine $1 without the bias of the 

264 global momentum conservation effect. One possibility 

265 might be to use the modified event plane method from 

266 Rcf. [20]. 

267 To get a feeling on how many V£* m terms are needed to 

268 exhaust the information encoded in distribution fc($* — 

269 $*„), in Fig. 3 we show the centrality distribution of V^* m 

270 for several values of j for j > 1 . Most of correlations are 

271 exhausted by including the j = 1-5, except for a few 

272 cases at N paxt < 100, such as 4($| - <f>|) and 3($| - $3). 

273 The results presented in Figs. 1-3 are obtained using 

274 the r 2 weighting (i.e. Eq. 2) and Glauber model. Altcr- 

275 natively we have calculated the $„ using the r n weighting 

276 for n > 1 and r 3 weighing for n = 1 [3] ; we also repeated 

277 the same calculation using a CGC (Color Glass Conden- 

278 sate) geometry [28] for both the r 2 and r n weighting. 

279 The results for these three cases are shown in Fig. 4 for 

280 j = 1. The main observations are qualitatively similar 
28i to Fig. 2. However, there are some interesting changes 

282 in the magnitudes of some correlation values: the use of 

283 CGC model increases the correlation for (71,777) =(2,4) 

284 and (2,6) but decreases the correlation for (71,777) =(3,6) 

285 in mid-central collisions; the use of r n weighting also in 

286 general increases Vj* m and affects the relative magnitude 

287 of Vj[* m between (77,771) =(2,4) and (3,6) in central colli- 



Due to the phase shift between the two types of cor- 
relations given by Eq. 25, many positive V^* m values in 
Fig. 2 and 3 would imply the corresponding V^ m values 
are negative. This happens for odd j and (n,m) =(2,3), 
(3,6), (2,5), (1,2), (1,4) and (1,6). If the n th -order flow di- 
rection align with $„ as predicted by the Glauber model, 
one should expect the signs of the correlators between the 
experimental event plane in Eq. 6 to exhibit very interest- 
ing dependence on choice of (71,771). On the other hand, 
if the dynamic mixing between flow of different orders 
is important [9], then this dependence could be strongly 
distorted. Therefore, direct measurements of the cor- 
relations between the experimentally determined event 
planes of different order can help to resolve this issue. 
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FIG. 3: (Color online) Centrality dependence of higher-order Fourier coefficients of the correlation Vn* m = (cos jfc(^ — ^m)) 
for various choices of n and m for j — 2 — 5 (from left column to right column). 
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FIG. 4: (Color online) Centrality dependence of first-order Fourier coefficients of the correlation V^* m = (cos k^^ — ^m)} for 
various choices of n and m for different ways of calculating e n (from left column to right column). 



B. Correlation between three planes 



304 It is straightforward to carry out the study of corre- 

305 lations between three planes. The "*" notation is again 

306 used to indicate the plane calculated using the major 

307 axes. The top panels of Fig. 5 show the 2-D normalized 

308 distribution d 2 N cvts / ^$2,1^3,1) in 40-50% centrality 

309 interval; the corresponding V^'l 3 coefficients are shown 

310 in the bottom panels. The coefficients along i = or 

311 j = simply reflect two plane correlators, V(* 2 and V^", 

312 respectively. The interesting coefficients are those for 



313 i,j ^ 0. A tight diagonal correlation in the top panels 

314 can be identified with a large (i,j) = (1, —1) component, 

315 which corresponds to (cos($ 1 + 2$2 — 34>g)}. This cor- 

316 relation is very weak in central collision and increases 

317 gradually towards peripheral collisions (see Figs. 11-14), 

318 similar to the finding in [3] (there is a sign difference 

319 due to the use of major axes here). This correlation is 

320 also observed to be generally bigger for r" weighting and 

321 CGC geometry. Hence a precise determination of this 

322 correlator could allow us to distinguish between different 

323 models for initial geometry. Sizable coefficients are also 

324 observed for (i,j) =(1,-2), (2,-2) and (1,1), correspond- 
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FIG. 5: (Color online) (top row) The correlation (normalized to 1) between $2,1 an d ^31 in 40-50% centrality interval for 
different weighting (r or r n ) and initial geometry models (Glauber or CGC). (bottom row) The corresponding Fourier spectrum 
1^*2*3. The constant term (i,j) = (0,0) (or V 1 ' 2 3 = 1) is omitted for clarify. The typical statistical error on the Fourier coefficients 
is about 0.002. 



-2$£-6$3)>, (cos(2$|+4$^-6$g)) and 
- 3^3)), respectively. Also note that the 
327 coefficients for (i,j) = (2,-1) and (3,-2), corresponding to 



325 ing to (cos(4$j 
(cos(5$t - 2$^ 



(cos($* 



4$^ 



3*3)) and (cos6($3 - $|)), are nearly 



329 zero consistent with the findings in Ref. [3] and Fig. 2, 

330 respectively. 

331 More results on other types of three plane corrcla- 

332 tions arc summarized in Figs. 6-10. It is generally 

333 observed that the Fourier components are always big- 

334 ger for r n — weighting than for r 2 —weighting, and they 

335 are slightly bigger for the CGC geometry than for the 

336 Glauber geometry. Some of the correlators are quite 

337 large, e.g. (cos(2$|+2$|-4$|)), (cos(2$t-6$^+4$|)), 

338 (cos(2$I+6$^-8$|)), (cos($t+3*3-4$|)}, (cos(2$J- 

339 6$|+4^)),(cos($J+4$|-5^)), (cos(3$t+2$2- 5$ 5))> 

340 (cos($t - 6$^ + 5$g)>, (cos(2$^ + 3$^ - 5$|)), and 

341 (cos(4$j — 9$3 + 5^5)). These correlators are shown as 

342 a function of centrality in Figs. 11-14. In general, they 

343 all increase from central to more peripheral collisions, 

344 however the rate of the change depends on the type of 

345 the correlator. The correlator (cos($5; + 2<f>2 - 3$g)) has 

346 the largest values in most cases, except for r" weighting 

347 in central and mid-central collision, where the correlator 

348 (cos(3>* + 3$3 - 4$|)) has the largest values. 

349 Interestingly, most of these correlators, when defined 

350 relative to the major axis, are positive (the negative val- 

351 ues are indicated with red "x" in Figs. 5-10). However, 

352 some of these correlations are likely to be strongly dis- 

353 torted due to the mixing during the hydrodynamic evolu- 

354 tion, especially for those involving $| and $5 . Neverthe- 

355 less, it would be interesting to measure these quantities 



experimentally and compare with our predictions. 



357 IV. DISCUSSION AND CONCLUSION 

358 In summary, we discussed a method for measuring the 

359 correlations between the directions of the anisotropic flow 

360 of different orders. This method involves Fourier trans- 
36i forming the differential distribution of the correlations 

362 between different event planes into various Fourier com- 

363 ponents, where each component is corrected separately 

364 by an event plane resolution term. This method has 

365 the advantage of simultaneously analyzing many diffcr- 

366 cnt correlators, especially those involving three or more 

367 different event planes, thus help identifying significant 

368 components. 

369 The expected strength of various two- or three-plane 

370 correlators are estimated using a Monte Carlo Glauber 

371 model or CGC model. Strong positive correlations are 

372 seen between the major axes of two eccentricities in 

373 mid-central and peripheral collisions for (n,m) =(2,4), 

374 (2,6), (3,6), and those involving the first-order eccentric- 

375 ity. Similarly, several significant three-plane correlators 

376 have also been identified, revealing novel correlation pat- 

377 terns expected from the average geometry and/or initial 

378 state fluctuations. These strong correlations imply the 

379 need to measure several two- or three-plane correlators 

380 in order to describe the full distribution. A detailed com- 
38i parison of the correlations calculated here with the data 

382 could shed light on the role of the initial geometry fluc- 

383 tuations and dynamic mixing during the hydrodynamic 
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FIG. 6: (Color online) The correlation between $2,1 an d $4,1 (top row) and the corresponding Fourier coefficients (bottom 
row) in 40-50% centrality interval. 
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FIG. 7: (Color online) The correlation between $y and $4_i (top row) and the corresponding Fourier coefficients (bottom 
row) in 40-50% centrality interval. 
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FIG. 8: (Color online) The correlation between $| ]2 and $4_2 (top row) and the corresponding Fourier coefficients (bottom 
row) in 40-50% centrality interval. 
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FIG. 9: (Color online) The correlation between $2,1 and $5,1 (top row) and the corresponding Fourier coefficients (bottom 
row) in 40-50% centrality interval. 
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FIG. 10: (Color online) The correlation between Sy and $5^ (top row) and the corresponding Fourier coefficients (bottom 
row) in 40-50% centrality interval. 



10 




r 2 Glauber^ K'^'^ 

» (cas(4&*+2<i> 2 *-G& 3 *)) 
o — >&— <cos(24> ] *+4* 2 *-&J> 3 ')>~ 

- (cos(5* 1 *-2* a *-3* 3 *)) 



fro-* 



- (cos(* 1 *+3*3*-4* 4 *)> 

- <cos(2* 1 *-6* 3 *+4* 4 *)> 

- (cos(* *+4* 2 *-5* 5 *)) 

- (cos(34> ] *+2* 2 i -5<I> 5 ')) 



i a 



. (cos(2* ] *+2* 2 *-4<V)> 
- (cos(2* ] *-6# 2 '+4 < V)> 
. (cos(2* *+6* 2 '-8<V)> - 



s ^ 



- (cos(<J> *- 

- (cos(2* ] 

- (cos (4* 



S# 3 *+5# 5 *)) 

+ 3* 3 '-5* 5 ')) 

-9<V+5* 5 ')> 



<&* o 



■=■ & 



-4 — *- 



Rl r 2 cgc "~ s>— < COS (*r +2 *2*- 3 *3*)> 

_^_ (cos(2<&*+4* 2 *-64> 3 *))- 
o 
o — B— (C0S(» 1 *-2® 8 *-»3 , )> 



_ □ * 



0.6 £ 



_S i_ 



1° on # 



- (cos{<£ *+3* 3 "-4<I> 4 *)> 

- <cos(2* 1 *-6* 3 *+4* 4 *)> 

- (cos(* *+4* 2 *-5* 5 *)) 

- (cos(3* 1 *+2* 2 *-5* 5 *)) 



ft 8 



. (cos(2* ] 
- (cos(2* ] 
. (cos(2* 



+2* 2 *-4<V)> 
-6* 2 "+44> 4 *)> 
+6* 2 '-8<V)> - 



# o 

vj& o 

1 5lo 



- (cos(*"-i 

- (008(2®,* 

- (005(4** 



* 3 *+5* 5 *)} 

h3* 3 *-5* 5 *)> 

9*3*+5*b*» 



ft «t 

e 4 



J — * & 6- 



FIG. 11: (Color online) The centrality dependence of vari- 
ous significant three plane correlators for r 2 weighting with 
Glauber geometry 



FIG. 13: (Color online) The centrality dependence of various 
significant three plane correlators for r 2 weighting with CGC 
geometry. 
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FIG. 12: (Color online) The centrality dependence of vari- 
ous significant three plane correlators for r n weighting with 
Glauber geometry. 



FIG. 14: (Color online) The centrality dependence of various 
significant three plane correlators for r™ weighting with CGC 
geometry. 



